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' In this paper we propose a scheme for quasi-perfect state transfer in a network 

^ ■ of dissipative harmonic oscillators. We consider ideal sender and receiver oscillators 



connected by a chain of nonideal transmitter oscillators coupled by nearest-neighbor 
resonances. From the algebraic properties of the dynamical quantities describing the 
evolution of the network state, we derive a criterion, fixing the coupling strengths 
between all the oscillators, apart from their natural frequencies, enabling perfect 



I state transfer in the particular case of ideal transmitter oscillators. Our criterion 



provides an easily manipulated formula enabling perfect state transfer in the special 
case where the network nonidealities are disregarded. By adjusting the common 
frequency of the sender and the receiver oscillators to be out of resonance with that 
of the transmitters, we demonstrate that the sender's state tunnels to the receiver 
oscillator by virtually exciting the nonideal transmitter chain. This virtual process 
makes negligible the decay rate associated with the transmitter line on the expenses 
of delaying the time interval for the state transfer process. Apart from our analytical 
results, numerical computations are presented to illustrate our protocol. 



PACS numbers: 42.50.-p; 42.50.Ex; 42.50.Lc 



I. INTRODUCTION 



A great deal of attention has been devoted recently to the subject of perfect state trans- 
fer (PST) in quantum networks. Since an actual quantum processor would require, in fact, 
the ability to transfer quantum information between spatially separated interacting systems 
composing a network, protocols have been established for PST in many different and gen- 
eral contexts. Among the several interesting theoretical contributions, aiming to advance the 
understanding and eventual implementation of a quantum processor, the construction of ef- 
fective two-qubit gates from PST between distant nodes in engineered bosonic and fermionic 

n 

networks was proposed p|. A general formalism of the problem of PST in networks of any 
topology and coupling configuration was also developed 2|]. Focusing on spin chains, PST 
has been pursued in networks extending beyond the nearest-neighbor couplings [sl, and a 
class of qubit networks allowing PST of any state in a fixed period of time has been devised 



1^. The problem of the scaling of errors (arising from network nonidealities) with the length 
of the channel connecting the nodes has been discussed j^. In this connection, a protocol 
for arbitrary PST in the presence of random fluctuations in the coupling strengths of a spin 
chain was reported j^. The entanglement dynamics in spin chains subject to noise and dis- 
order has also been analyzed in Ref. [3]. Regarding PST in networks of harmonic oscillators, 
a comprehensive analysis of this subject has been presented in Ref. j^, and a protocol for 
high-efficiency transfer of quantum entanglements in translation-invariant quantum chains 
has been proposed 9|]. In contrast to the achievements with spin chains, there have been 
no proposals, until this paper, for arbitrary PST through a network of nonideal harmonic 
oscillators. 

Since PST is achieved by appropriate tuning of the intermode coupling and the frequencies 
of the systems composing the network, the search for the general rules governing these 
adjustments is of crucial interest. In this regard, it is worth noting that the general recipe 
presented in Ref. ^ contrasts with the procedure in Ref. 9| where a degenerate chain of 
oscillators is considered: all oscillators having the same frequency and interacting with each 
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other with the same couphng strength. The lack of PST in Ref. is compensated by a 
more reahstic possibihty of implementation of a translation-invariant network, apart from 
the high efficiency transfer operation reported. We observe that in the present contribution 
we also derive general conditions for PST, with the advantage that they are easier to handle 
mathematically than those outlined in Ref. [2]. 

As well as in spin chains, significant advances have been made recently in bosonic networks 
A general treatment of a network of coupled dissipative quantum harmonic 
oscillators has been presented recently , for an arbitrary topology, i.e., irrespective of the 
way the oscillators are coupled together, the strength of their couplings, and their natural 
frequencies [l^. Regarding the dissipative mechanism, two different scenarios are considered 
in Ref. [l^: a more realistic in which each oscillator is coupled to its own reservoir and 
another with all the network oscillators coupled to a common reservoir. Within such a 
general treatment of dissipation, the emergence of relaxation- and decoherence-free subspaces 
in networks of weakly and strongly coupled resonators has also been addressed 15|]- We 
finally mention the proposition of a quantum memory for the preservation of superposition 
states against decoherence by their evolution in appropriate topologies of such dissipative 
bosonic networks [l^. 

Building on the achievements mentioned above, in the present study we develop a protocol 
for quasi-perfect state transfer (QPST) in a linear network of dissipative oscillators. More 
specifically, we envisage the transfer of a state between ideal sender and receiver oscillators 

nn 

through a linear chain of nonideal transmitter oscillators. By analogy with Refs. p^, I15| . 
here we again adopt the general and more realistic scenario where each transmitter oscillator 
is coupled to a distinct reservoir. Therefore, our goal resembles that of arbitrary PST in 
Ref. gI; however, instead of fiuctuations in the couplings of a spin chain, we deal with 
the fiuctuations injected by each of the reservoirs coupled to the transmitter oscillators. 
Anticipating our strategy to achieve QPST despite these sources of nonideality, we adjust the 
frequency of the sender and the receiver oscillators to be significantly out of resonance with 
that of the transmitters; within such an arrangement the state to be transferred occupies 
the transmitter oscillators only virtually, weakening the undesired effects of their decay 
mechanisms. 

We point out that the recent proposal of a variety of resources for wiring up quantum 



systems 17|] lends a strongly realistic bias to the possibility of controlling the transfer of 



information in quantum networks. Following the mastering of the manipulation of the in- 
teraction between single atoms and vibrating modes of high-Q cavities and trapped 
ions [l^ , circuit cavity quantum electrodynamics and photonic crystals seem to enhance the 
ability to transfer quantum information to a level enabling the implementation of a logic pro- 
cessor {2^. In this connection, the elaboration of schemes to circumvent the noise injection 
in the processes of PST in networks of coupled nonideal quantum systems is indispensable, 
enabling protocols for fault-tolerant information transfer and deepenning our understanding 
of fundamental quantum phenomena such as entanglement and decoherence. 

The plan of the paper is as follows. In Section [Tll by analogy with the developments 
presented in Ref we introduce our model in which two ideal oscillators are connected 
by a transmission line of nonideal oscillators. The master equation describing the network 
dynamics is presented, together with its solution. In Section IIIH we present a criterion for 
PST (QPST) in ideal (nonideal) networks of harmonic oscillators, together with a particular 
application where the sets of parameters {iVm} and {Xmn} ensuring PST are derived. Our 
criterion relies on the definition of matrix 0(t), describing the evolution of the system, which 
serves two purposes: to ensure PST and to compute the exchange time tcx- The analytical 
treatment of QPST in small nonideal linear networks of 3, 4, 5, and 6 oscillators are given 
in Section ??. In section |V] we apply numerical procedures to extend our results to large 
numbers of nonideal transmitters, providing a comprehensive analysis of all the network 
parameters involved and demonstrating the robustness of our model. Finally, we present 
our concluding remarks in Section IVIi 



II. THE MODEL, THE CORRESPONDING MASTER EQUATION, AND ITS 

SOLUTION 

A. General bosonic dissipative network 

Before introducing our model for PST, we first revisit the developments in Ref. jl^ . 
considering a network of interacting dissipative oscillators from the general perspective 
where each oscillator interacts with each other. Any particular network topology (or graph) 
follows from this general approach, with an appropriate choice of the parameters defining 
the Hamiltonian modelling the network. As stressed above, by topology is meant i) which 



resonators are coupled together, ii) their couphng strengths and in) their natural frequencies. 
Moreover, in a more realistic approach for most physical systems, it is assumed that each 
network oscillator interacts with its own reservoir, instead of the special case where all the 
oscillators interact with a common reservoir. Therefore, assuming from here on that the 
indexes m, n, and i, labeling the oscillators, run from 1 to N, we start from the general 
Hamiltonian 



(1) 



where bl^f^ {hmk) is the creation (annihilation) operator for the fcth bath mode Umk coupled 
to the mth network oscillator u;^, whose creation (annihilation) operator is {dm)- The 
coupling strengths between the oscillators are given by the set {Xmn}-, while those between 
the oscillators and their reservoirs by {Krtfe}. The reservoirs are modeled as a set oi k = 
1, . . . , oo modes, and the elements Hmn defining the network topology compose the matrix 
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To obtain the master equation of the network we first diagonalize the Hamiltonian H 
through a canonical transformation Am = C*mn«n, where the coefficients of the mth line 
of matrix C define the eigenvectors associated with the eigenvalues of matrix ([2]). The 
commutation relations [^1^,^^ = Smn and = 0, following from the orthogonality 

of matrix C, in that = C~^, enable us to rewrite Hamiltonian ([1]) in terms of decoupled 
normal-mode oscillators Wm, each of them interacting, however, with all the reservoirs. 
The diagonalized H helps us to introduce the interaction picture, within which a set of 
assumptions leads to the master equation describing the network evolution. We first assume 
that the couplings between the normal-mode oscillators and the reservoirs are weak enough 
to allow a second-order perturbation approximation. Moreover, under Markovian white 
noise, the evolved reduced density operator of the network, p{t), is assumed to be factorized 
from the stationary density operator of the reservoirs. Finally, the reservoir frequencies are 



(2) 



assumed to be closely spaced enough to allow a continuum summation, the spectral density 
(Tm(w„) and coupling parameter Vm{zUn) being slowly varying functions. Thus, after tracing 
out the degrees of freedom of the absolute zero reservoirs, we obtain the generalized Lindblad 
form 

m,n ^ 
P 

+ ^ {[anp{t),al^] + [am,pit)al]) 

= 5Z ) ^ [P(^), al^H^nan] + Cmnp{t) > , (3) 
m,n ^ ■' 

where we have defined the effective damping matrix F, whose elements are 

n' 

with 7m (^n) = Jf [Vm{ )] 5{e)de. In Eq. ([3]), the Liouville operators 

^mnP{t) account for both the direct (m = n) and indirect (m 7^ n) dissipative channels. 
Through the direct dissipative channels, the oscillators lose excitation to their own reser- 
voirs, at a damping rate Tmm, whereas through the indirect channels they lose excitation 
to all the other reservoirs but not to their own. For Markovian white noise reservoirs the 
indirect channels disappear, since the spectral densities of the reservoirs are invariant over 
Tanslation in frequency space, rendering 7m(w„') = 7^, and, consequently, = ^7m<^mn 
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22|. 



It is worth stressing that the whole of the working presented in the next subsection, where 
a solution to the master equation is derived, applies to the case of non-Markovian reservoirs, 
thus including the indirect dissipative channels. Markovian reservoirs are only assumed at 
the end of subsection C, when the topology of our framework is finally defined as a linear 
chain of dissipative oscillators. 



B. Solution of the master equation 

To obtain a solution of Eq. ([3]), we shall employ it to derive the Glauber-Sudarshan P 
function for the network 



where we have defined the matrix elements = iHmn + rmn./2, extending the former 
Hmni to take into account the dissipative {D) process. Assuming that the general network 
described by matrix ([2]) is composed entirely of dissipative oscillators, the matrix assumes 
the form 

Til ri2 • 

r22 ■ 
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For the initial state of the network we consider the general pure superpositions of coherent 
states 

Q 

p{Q) = KK\Wm}) . (6) 

r,s=l 

where M is the normalization factor, is the probability amplitude of the product state 
|{/^m}) — ®m=i l/^m)' ^^'^ the labels r and s run from 1 to the integer Q. The superscript r 
stands for the rth state of the superposition while the subscript m stands for the coherent 
state of the mth oscillator. (We stress that the discrete sum of product states in Eq. ([6]) 
can be substituted by the continuum sum |'?A(0)) = Af f d9A {6) \ {[3m {0)}) with no further 
complication). From Eqs. (jll) and (I6l), it is straightforward to show that the network density 
operator evolves as 



(7) 



({awiKcw}) 

The excitation of the mth oscillator, given by 

Cmit) = Y.^mn{t)[^l (8) 

n 

follows from the time-dependent matrix elements 

0mn(t) = ^ Drnm' exp {-Wm't) D~}^. 

m' 

where the mth column of matrix D defines the mth eigenvector associated with the eigen- 
value Wm of matrix H^. 

For the reduced density operator of the mth oscillator we obtain 

(9) 



where the infiuence of all the other oscillators of the network is present explicitly in the 
product {{Pn} \{Pn}) and implicitly in the states |Cm(^))- 



C. A linear dissipative network: our framework 



A linear (lin) dissipative network of nearest-neighbor interacting harmonic oscillator is 
built up by coupling the kth oscillator with the {k ± l)th oscillators, leaving the first oscil- 
lator (m = 1) uncoupled from the last one (m = A^). The matrix H^^ obtained for this case 
has the three-diagonal form 
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As already made clear, we focus on the case of ideal sender and receiver oscillators, here 
assumed to be the first and the last, respectively, both with the same frequency u. All 
the transmitter oscillators, from the second to the {N — l)th, are assumed from here on 
to decay at the same rate Tmn = T, and to be tuned out of resonance with the sender 
and receiver, to frequency Q. Regarding the coupling between the oscillators, we assume 
that the sender and the receiver are connected with their transmitter neighbors with the 
same strength A, whereas the transmitters are connected to each other with strength eX, e 
being a dimensionless parameter allowing the couplings within the transmitting channel to 
be controlled. Finally, assuming Markovian white noise reservoirs, to eliminate the indirect 
dissipative channels, the above matrix simplifies to 

(uj \ ■■■ 

A VL-iVl2 e\ ■■■ 
e\ VL-iV/2 ■•. 
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III. A CRITERION FOR PST (QPST) IN IDEAL (NONIDEAL) NETWORKS 

OF HARMONIC OSCILLATORS 



A. PST in ideal networks 

On the basis of the model developed above, in this section we derive a general criterion 
for PST, whatever the topology of an ideal (F = 0) network of harmonic oscillators. Such 
a criterion for PST is further applied for the particular case of the linear topology modeled 
by Eq. ( ITOl) . Starting from Eq. ([8]) written as the matrix product (t) = @(t) ■ (3^ , it is 
straightforward to conclude that the condition for transferring the state of the first oscillator 
to the A^th is given by the matrix structure 

/ n n ... i\ 
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(11) 



which reversely ensures that the state of the A^th oscillator can equally be transferred to 
the first one, no matter what happens to the states of the intervening N — 2 transmitting 
oscillators, from the second to the last but one. That is why the submatrix formed by 
removing the first and last rows and columns of is left undefined; for any submatrix [V], 
PST is achieved as long as the states of the first and the last oscillators are interchanged. 
It is worth noting that the criterion for state exchange fixed by the structure of the matrix 
in Eq. (fTTj) does not depend on the network state. Evidently, the imposition of this matrix 
form can be used to compute the exchange time tex- 
Since the evolution of matrix 0(t), given by 



&{t) = D-exp[-zRt] ■ 

= exp[-(D-iR ■ D~^)t] = exp[-iHt], 



(12) 



is governed by the Hamiltonian H (whose eigenvalues Rm compose the diagonal matrix R), 
the sets of parameters {um} and {Xmn} complying with a necessary condition for PST in a 
particular network topology follow from the commutation relation 



[0(tex),H]=O, 



(13) 



which holds at any time, including the state exchange time tex- We observe that not all sets 
{uJm} and {Xmn} derived from condition ffT^ ensure PST. For a necessary and sufficient 
condition we must choose, among the sets of parameters {ujm} and {Xmn}, those ensuring 
the reduction of matrix ( fT2l) to ( fTTj) . our initial premise, at the exchange time. Hence, matrix 
@{t) can be used for two purposes: to ensure PST and to compute the exchange time tex- 

We observed in the Introduction that our general condition for PST, given by Eq. (IT^l. 
has the advantage of being easier to handle mathematically than that outlined in Ref . [2 1 . 
In fact, in Ref. 2|], the condition for PST follows from the diagonalization of a permuta- 
tion matrix — equivalent to the above defined 0(t) — whereas our condition reduces to the 
computation of a commutation relation. 

It must be stressed that in the case where the state transfer is achieved by tunneling, the 
set of coupling parameters {Xmn} does not play a significant role in the process since the 
transferred state does not effectively occupy any of the transmitter oscillators. Actually, for 
the state to be transferred by tunneling, we only have to ensure that the common frequency 
(uj) of the sender and the receiver oscillators is distinct from those of the transmitter ones; 
under this condition, we can verify the relation ( fT3l) independently of the set of coupling 
strengths Xmn, and these may all have the same value A. 



B. QPST in nonideal networks 

The extension of the above criterion to state transfer in nonideal networks, the case at 
hand, follows directly by assuming a finite value of P. The damping rates due to dissipation 
are thus introduced into Eqs. (ITT]) and (|T^ . generalizing the relation (|T^ to 

[e(tex),H^] =0. (14) 

However, in this general nonideal case the transferred state exhibits a fidelity that decreases 
with time and the magnitude of the decay rate P, apart from depending on other properties 
such as the network topology. Moreover, there are two further situations, not envisaged in 
this paper, where the commutation relation ( IT^ reduces to that of the ideal case (|T3i) . even 
when the fidelity decreases with time and the magnitude of P: i) when Tmn = Tf^mn and ii) 
Tmn = T- In both cases the sets of parameters {tOm} and {Xmn} ensuring QPST are exactly 
those ensuring PST. 



It is worth noting again that, confining ourselves to a hnear dissipative network, the 
adopted strategy of considering the ideal sender and receiver oscillators significantly out 
of resonance with the nonideal transmitters, ensures QPST despite the nonidealities of the 
transmitter line. In fact, as anticipated in the Introduction and demonstrated below, the 
virtual occupation of the transmitter line protects the transferred state almost perfectly 
from the damping mechanisms. 



C. An application of the commutation relation [0(tex),H] = 



Considering the Hamiltonian H in Eq. ([2]) and the particular choice 

/O ■•■ l\ 



(15) 
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1 ■■■ 
\1 ■•■ 0/ 

ensuring the transfer of the initial state of the mth oscillator to the [N — [m — l)]th one, 
we obtain from relation ffT3l) the condition for PST: 



— ^N-{m-l)^ 
Amn. = ^N~{m.~l),N-{n-l), 



which generalizes the one introduced in Ref. |8|, given by 
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— m). 



(16a) 
(16b) 



Evidently, each choice of the submatrix [V] in the general form ffTTj) prompts a different set 
of parameters {uJm} and {Xmn} ensuring PST. 



D. An alternative way to compute tex 

Another way to compute the exchange time tex? instead of using matrix G(tex); follows 
from the time profile of the probability of a successful transfer of the desired state — or 



equivalently, the fidelity of tlie transfer process — given by 

nt) = Tr [p(0)p(t)] . (17) 

Evidently, in a general situation where PST is sought — other than the transfer by tunneling 
analyzed here — the computation of the exchange time tox from Eq. (1171) . must take into 
account the set of parameters {oOm} and {A^n} ensuring PST, derived from the commutation 
relation ([Hj). A further condition is that the ideal case F = must render 'P(tcx) = 1- 
Although a rather demanding task, the maximization of the probability 'P(tcx) yields, by 
itself, a necessary and sufficient condition for the derivation of the sets of parameters {oOm.} 
and {Xmn} ensuring PST. The commutation relation (I14p provides a shortcut for this task. 

We assume, as usual, that the initial density operator p(0) factorizes as Pi{0)l^^_2Pm{0) , 
pi(0) representing the state of the first oscillator, to be transferred to the A^th one, and 
Pm(0) standing for the initial state of the transmitter and the receiver oscillators. In its 
turn, the density operator at the exchange time p(tox) factorizes as P7v(tcx)(S)m=iPm(^ex), 
where pAr(tcx) represents the transferred state of the A^th oscillator, and Pm(^cx) the final 
state of the transmitter oscillators plus the sender. Evidently, the state of the A^th oscillator 
at tcx must be, in the ideal case where F = 0, exactly that prepared in the first oscillator. 
With this assumption, by substituting Eq. ([9]) into 'P(tcx), we obtain 



/3, 



X exp 



cx) Pl\ Ca''(^cx) Pi 



(18) 



As an illustrative application of Eq. ( fTSi) . we consider the specific case where the state 
to be transferred to the A^th oscillator is prepared in the first one as the Schrodinger cat- 
like superposition M {\a)^ + \—a)^), while all other oscillators are in the vacuum state. We 
obtain from Eq. ffTSl) . for the ideal case, the relations 



'^Cntyi. cos {Rmtcx) — ±1, 

m 

Yli^Nm sin {Rmtex) C'ml = 0) 



(19a) 
(19b) 

m 

which enable us to determine tex as well as Eq. ( |T3l) . remembering that Rm stands for the 
eigenvalues of the free Hamiltonian H. Therefore, the exchange time tgx follows from the 
computation of the eigenstates and the associated eigenvectors of H. 



IV. ANALYTICAL TREATMENT OF QPST IN SMALL NONIDEAL LINEAR 

NETWORKS 



Now, focusing on our scheme for the transfer of states by tunnehng , we present an 
analytical treatment of small nonideal linear networks, using the particular case of ( ITSl) to 
compute the exchange time t^x- As an illustration, consider the case of QPST in a network 
with = 4, whose dissipative matrix, following from Eq. (IIUI) . is 

fiuj iX 0\ 

iX in + T/2 zeX 
ieX in + T/2 iX 
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\0 iX iujj 

The eigenvalues and eigenvectors of H^^j.^^^^ defines the matrix 0(t) and, consequently, the 
evolved excitation of the oscillators C^' (t) = &{t)-(3^ . Next, we introduce the scaled damping 
rate 77 = F/A, frequency w = uj /X, and time r = At, apart from the detunings A-t = f2 ± u; 
and the effective coupling = A/A_, to stress that we must focus on the regime where 
/i ^ 1 and ?7 ^ 1. These parameters ensure the weak coupling between the transmitter 
oscillators and their respective reservoirs — ^justifying the master equation derived above — 
apart from enabling the expansion of matrix 0(t) to second order in /i, giving 

/ cos(£:;uV) -Qc{r) iQs{r) -z sin(£:/iV)^ 

-Bc(r) f)c(r) -if).(r) is,(r) 

iQs{,r) -i^s{,r) f)c(r) -0c(r) 

\-2sin(£:/iV) 2g,(r) -Qc{.t) cos(£:/iV) J 

where we have defined the functions 
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o{^?), (20) 



0c(t) = /i [cos(£:/iV) - fic(r) 
Qs{t) = jj, [sin(£:/iV) - [),(r) 

= exp [- [t] + r] cos (er) 
f)s(r) = exp [- [t] + z/i"^) r] sin[(£:r) 



From the condition for transfer of the state of the first oscillator to the last one, fixed 
by the matrix structure (|T5l) . and assuming the additional restriction (efi)'^ <^ 1, we thus 



obtain from Eq. (1201) the scaled exchange time 



which imphes the relations 



CI (rif =^)) ^ e-^'^-- [PI - {^Qc{rj:^W2 + Qsir^^m)] + O (^^) , (21a) 
CI (Tir^)) - e-^-^-'^^ [(31 - {Qs{r!:^W, + ^Qc{r!:m)] + O (^^) , (21b) 

to ensure that the commutation [0(rex ),H^] = is satisfied. 

As is evident from the above expressions, the relaxation process represented by the scaled 
damping rate rj prohibits a perfect state transfer, attenuating the value of the excitations 
PI and PI- Moreover, the reservoirs also spoil the desired relation Q^^-^ ^rijf"^'' j = Plf^^^ by 
mixing it with the excitations of the nonideal transmitter oscillators. However, and this is the 
core of our technique, the tunneling mechanism of state transfer prompts the decay function 
g-ry/i^rix' _ g-7rr/2eA approachcs Unity — within the ranges of the parameters outlined 

above, i.e., /i, 77, (e/i)^ <^ 1 — as the couphng eX between the transmitting oscillators is 
increased. We note that the increase of eX decreases the exchange time Tex~^'^ — 7rA_/2yU£:A, 
and that is the reason for the choice of a strong coupling strength sX between the transmitter 
oscillators: to enable the control of the exchange time Tex, decreasing it as the fidelity of the 
transfer process increases. Otherwise, without the dimensionless parameter e, increasing of 
the detuning A_ (i.e., the process allowing a significant fidelity of the transferred state), 
would result in an uncontrolled rise in the exchange time Tex~^^ — 7rA_/2/iA. From Eq. 
([21]), we conclude that the fidelity of the state transfer mechanism is maximized when the 
transmitter oscillators are prepared in the vacuum state, apart from weakening the inevitable 
system-reservoir coupling. In fact, the amount of excitation of the dissipative transmission 
channel is directly proportional to the intensity of noise injected into the transferred state. 
We finally note that for the case of an ideal transmission channel, i.e., 77 = 0, we obtain — 
up to first order corrections in carried out in gcl^cx^) and Qs{tcx) in the general case of an 
excited transmission channel — the desired relations Q ^ri^^"^'' j ^ PI and Q ^ri^f j — 

On the basis of the exchange time Tex = 7r/2 for the simplest network, composed of two 



oscillators with coupling strength A 22[, it is useful to assign an effective coupling strength 



to the network which, for the case = 4, turns out to be A^^ ~ efi'^X. 



A. Nonideal linear network of = 3, 5, and 6 oscillators 



Following the steps described above for the case = 4 and adopting the same regime of 
parameters fj,,ri, (e/i)^ ^ 1, we obtain for = 3 the result 

giving the expected effective coupling strength A^^"^"* ~ 2fiX between the sender and re- 
ceiver oscillators. Evidently, it is also possible to derive analytical expressions for the scaled 
exchange time for networks with > 4 whenever the diagonalization of the associated ma- 
trix H^^j.^^ generates a characteristic polynomial that factorizes into parts of degree < 4. 
We have found that = 8 is the limiting case permitting analytical solution, with a charac- 
teristic polynomial that factorizes into two parts of degree 4. For A^ = 9, the characteristic 
polynomial factorizes into two parts of degrees 4x5. Analyzing the cases A^ = 5 and A^ = 6, 
which factorize into polynomial of degrees 2x3 and 3x3, respectively, we obtain the results 



[1 + OQ/)] 



We finally note that, without the imposition of the restriction (e/i)^ <^ 1, the above 
expressions for Tcx would be approximately rewritten, for small values of A^, as 

{1+[A+B7]'-Ce']fi^)}, (22) 



' ex 



2e(^-3)^(iv-2) 

where ^ = A^ — 1, i3 = X]m=i^5 ^ '~ 3- Evidently, for small A^ and (^yu)^ <^ 1 the 
above expression is equivalent to the results derived from A^ = 3 to A^ = 6. The derivation 
of an analytical expression for Tc^^ in the general case of any A^ is not an easy task. As 
A^ increases, we find that an involved dependence of the second-order correction 0{fi^) on 
A^ begins to play a significant role. Although the task of identifying this dependence is 
still a compelling challenge, in the present paper we analyze QPST for large values of A^ 
numerically, using expression ( fTTl) to obtain the fidelity of the transfer process. 



V. NUMERICAL TREATMENT OF QPST IN LARGE NONIDEAL LINEAR 

NETWORKS 



Now, using the probability of a successful transfer of the desired state, given by Eq. (ITTj) . 
we analyze the QPST for large nonideal linear networks. Consider the transfer of state 
A/'(|a)-|^ + |— to the A^th oscillator, with a = 5 and all other oscillators in the vacuum 
state. In Fig. l(a, b, c, and d) we plot the numerical curves for the exchange probability 
Vcx{t) against r for the cases = 5,10,50, and 100, respectively. Considering the above 
regime of parameters, /i, ?7 ^ 1 — without requiring the additional restriction (efi)'^ ^ 1 — 
we assume, in units of the coupling strength A, the fictitious value w = u = 10, giving 
A_ = = 10^ and e = 5 x 10^, apart from rj = T = 10~^. We first observe, as expected, 
that the exchange time Tc^ '^ increases proportionally to A^, being around Te^'^ ^ vr x 10^ for 
the case of iV = 5, in agreement with the analytical result computed in Eq. ([22]). Moreover, 
as the state to be transferred occupies the virtual nonideal channel for a time interval 
proportional to N, the fidelity of the transfer process decreases with A^, as displayed in Fig. 
1. In fact, for the choice of decay rate outlined above, we show that the fidelity of the 
transfer process is about unity for the cases A^ = 5 and 10, beginning to exhibit a significant 
decrease from A^ = 50. We have not found any sensible decrease of the fidelity for the cases 
N = 5 and 10, even for time intervals many orders of magnitude longer than the exchange 
time. The shaded regions in the figures follow from the strong oscillations of the probability 
VcxiT), coming from the natural frequencies of the oscillators. 

Focusing on the case A^ = 10, in Fig. 2(a) we plot the exchange probability 'Pcx(t) 
against r, considering the same parameters as in Fig. 1, except for the excitation /5 = 5 of 
the coherent states assumed to populate the oscillators composing the transmitter channel. 
As expected, we verify that the exchange times are exactly those of the associated case in Fig. 
1(b). However, the pattern of the curves changes, giving rise to the expected background 
oscillations due to the initial excitation of the transmitter oscillators. Moreover, additional 
peaks occurs at the recurrence time Tree, when the superposition state Af{\a)-^^ + l~<^)i) 
goes back to the sender oscillator. The magnitudes of these secondary peaks follow from 
the probability of finding the superposition state in the receiver oscillator at the recurrence 
time, given by 'P(ri.cc) = Tr [pi(0)pAr(ri.oc)] = 1/2, where Pi(0) and Pn{ttcc) represent the 
superposition states standing in the first oscillator at t = and in the A^th one at tree- 



In Fig. 2(b) we assume the same parameters as in Fig. 2(a) to plot the probabihty of 
recurrence of the superposition state back to the first oscillator, given by the expression 
'Prcc(T) = Tr [pi(0)pi(r)]. As expected, in the recurrence time the superposition state recurs 
back to the sender oscillator. 

In Fig. 3(a) we adopt the same parameters as in Fig. 1(b), except for the smaller value 
e = 8 X 10^, to illustrate the expected increase in the exchange time. An interesting feature 
of this figure is that the fidelity of the transfer process does not depends on e, despite the 
increase in the exchange time, which is around 7 orders of magnitude greater than in Fig. 
1(b). These features are reinforced in Fig. 3(b) where, again with the same parameters as 
in Fig. 1(b), we consider the limiting case e = 1, which results in a formidable increase in 
the exchange time of around 28 orders of magnitude, still preserving the fidelity. 

With the same parameters as in Fig. 3(a), except for the smaller value of A_ = 2 x 10^, 
in Fig. 4(a) we first see that the smaller detuning pulls the exchange time back, to about 
the order of magnitude found in Fig. 1(b), in spite of the quantity e = 8 x 10^. Moreover, 
we verify the expected continuous decrease 

of the process due to such a small detuning, which compels the state to populate more 
effectively the virtual transmitter channel. In Fig. 4(b) we use the same parameters as in 
Fig. 4(a), except for the coherent state a = 10, to show that a larger excitation of the state 
to be transferred results in a smaller fidelity of the process. In fact, the decoherence time of 
a quantum state varies inversely with its excitation, in accordance with the correspondence 
principle. 

Finally, to illustrate the advantage of our tunneling-based scheme for state transfer over 
those where the transfer proceeds non-virtually through all the transmitter oscillators, in 
Fig. 5 we plot the exchange probability 'Pcx(t) against r for the case where the sets of 
parameters {ujm} and {Xmn} in Eq- dlSD are utilized. Assuming, as in Fig. 1(a), u = 10 and 
?7 = r = 10"^, we find that, whereas the fidelity of our tunneling scheme is about unity up 
to at least 10 times the exchange time TeK~^^ ^ 7 x 10^, that for the case of Fig. 5 decays 
to around zero for r ^ 2 x 10'^. 



VI. CONCLUDING REMARKS 



In this work we have discussed the problem of state transfer in a hnear chain of quantum 
dissipative harmonic oscillators. Assuming the first and last oscillators to be on-resonant 
with each other and significantly off-resonant with the transmitter channel — running from 
the second to the last-but-one oscillator — we take advantage of the tunneling effect to cir- 
cumvent the decoherence during the transmission process. We have assumed ideal sender 
and receiver oscillators connected by nonideal transmitters. Apart from the significant im- 
provement of the fidelity of our quasi-perfect transfer process compared to those where the 
transfer proceeds through a non-virtual mechanism, we have presented additional contribu- 
tions to state transferring processes. 

On the basis of a recently presented general treatment of a network of coupled dissipative 
quantum harmonic oscillators, we have also derived a general criterion for PST, whatever the 
topology of the network, which is applied here to the particular case of a linear chain. Our 
criterion, given essentially by the commutation relation between the network Hamiltonian 
H-^ and the matrix G(rcx) connecting the initial states of the network to those at the 
exchange time (when the transfer process is supposed to be accomplished) , has the advantage 
of being easier to handle mathematically than that outlined in Ref. In fact, the matrix 
0(t), defined in Ref. can be used for two purposes: to ensure PST (or QPST) and to 
compute the exchange time Tcx, which can also be controlled by manipulating the coupling 
between the transmitter oscillators. As a matter of fact, a large detuning between the on- 
resonant and the off-resonant oscillators ensures not only a high fidelity for the transfer 
process but also a significant prolongation of the exchange time. In this connection, an 
increase in the coupling between the transmitter oscillators acts to shorten such delays. 
The role played by each of the network parameters is analyzed in details in a set of figures 
presented after the formal development of our scheme. 

We thus propose that the tunneling scheme presented here can be useful for transferring 
quantum states between distant nodes of a quantum circuit without their undergoing a 
significant coherence decay. Evidently, the control of the network parameters, such as the 
natural frequencies and coupling strengths of the oscillators, still represents a sensitive issue 
to be overcome experimentally. We also point out that recent developments in circuit QED 
and photonic crystals signal realistic platforms for the experimental implementation of our 



proposed scheme. 

Acknowledgements 

The authors acknowledge the support from FAPESP and CNPQ, Brazihan agencies, and 
thank F. L. Semiao and G. D. de Moraes Neto for helpful discussions. 



[1] M.-H. Yung and S. Bose, Phys. Rev. A 71, 032310 (2005). 

[2] V. Kost'ak, G. M. Nikolopoulos, and I. Jex, Phys. Rev. A 75, 042319 (2007). 

[3] Alastair Kay, Phys. Rev. A 73, 032306 (2006). 

[4] M. Christandl, N. Datta, A. Ekert, and A. J. Landahl, Phys. Rev. Lett. 92, 187902 (2004). 
[5] H.-J. Briegel, W. Diir, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 81, 5932 (1998). 
[6] D. Burgarth and S. Bose, New J. Phys. 7, 135 (2005). 

[7] D. I. Tsomokos, M. J. Hartmann, S.F. Huelga, and M.B. Plenio, New J. Phys. 9, 79 (2007). 

[8] M. B. Plenio, J. Hartley and J. Eisert, New J. Phys. 6, 36 (2004). 

[9] M. B. Plenio and F. L. Semiao, New J. Phys. 7, 73 (2005). 
[10] K. Audenaert, J. Eisert, M. B. Plenio, and R. F. Werner, Phys. Rev. A 66, 042327 (2002). 
[11] A. S. M. de Castro , J. F. Peruzzo, and V. V. Dodonov, Phys. Rev. A 71, 032319 (2005). 
[12] M. Cramer, J. Eisert, M. B. Plenio, and J. Dreifiig, Phys. Rev. A 73, 012309 (2006). 
[13] C.-W. Chou, J. Laurat, H. Deng, K. S. Choi, H. de Riedmatten, D. Felinto, H. J. Kimble, 

Science 316, 1316 (2007). 
[14] M. A. de Ponte, S.S. Mizrahi and M.H.Y Moussa, Phys. Rev. A 76, 032101 (2007). 
[15] M. A. de Ponte, S.S. Mizrahi and M.H.Y Moussa, Annals of Physics (NY) 322, 2077 (2007). 
[16] M. A. de Ponte, S.S. Mizrahi and M.H.Y Moussa, J. Phys. B-Atom. Mol. Opt. Phys. 41, 
215506 (2008). 

[17] A. Blais, R-S Huang, A. Walhaff, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. A 69, 062320 
(2004); R. J. Schoelkopf and S. M. Girvin, Nature 451, 664 (2008); F. Illuminati, Nature Phys. 
2, 803 (2006). 

[18] J. M. Raimond, M. Prune, and S. Haroche, Rev. Mod. Phys. 73, 565 (2001). 

[19] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. Mod. Phys. 75, 281 (2003). 



[20] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cam- 
bridge University Press, Cambridge, 2000; The Physics of Quantum Information, edited by 
D. Bouwmeester, A. Ekert, and A. Zeihnger Springer, Berhn, 2000. 

[21] M. A. de Ponte, M. C. de Ohveira and M.H.Y Moussa, Phys. Rev. A 70, 022324 (2004); ibid. 
Phys. Rev. A 70, 022325 (2004). 

[22] M. A. de Ponte, M. C. de Ohveira and M.H.Y Moussa, Ann. Phys. (N. Y.) 317, 72 (2005). 
Figure captions 

Fig. 1. Numerical plots of the exchange probability 'Pcx(^) against r for the transfer of 
the state {\a)^ + \—a)^) to the A^th oscillator, with a = 5 and all other oscillators in the 
vacuum state. In units of the coupling strength A we set the fictitious value w = u = 10, 
prompting A_ = = 10^ and e = 5 x 10^, apart from rj = T = 10"'^. Curves (a)-(d) refer 
to the cases = 5, 10, 50, and 100, respectively. 

Fig. 2. (a) Numerical plots of the exchange probability Pcx(t) against r for the transfer 
of the state ^/ {\a)-^ + \ —a)^) to the A^th oscillator, assuming the same parameters as in Fig. 
1, except for the excitation /9 = 5 of the coherent states populating the oscillators of the 
transmitter channel, (b) Probability of recurrence VrcciT) of the initial superposition back 
to the first oscillator, plotted against r. 

Fig. 3. Numerical plots of the exchange probability Vcx{'t) against r for the transfer of 
the state M {\a)^ + |— a)-^) to the A^th oscillator, assuming the same parameters as in Fig. 
1(b), except for the smaller values of the parameter (a) e = 8 x 10^ and (b) e = 1. 

Fig. 4. (a) Numerical plot of the exchange probability 'Pcx(t) against r for the transfer of 
the state J\f {\a)^ + \ —a)-^ to the A^th oscillator, assuming the same parameters as in Fig. 
3(a), except for the smaller value of the parameter A_ = 2 x 10'^. (b) The same as in (a), 
but with the coherent state a = 10, to illustrate that the fidelity of the process decreases as 
the excitation of the state to be transferred increases. 

Fig. 5. Plot of the exchange probability 'Pcx(t) against r for the case where the network 
parameters in Eq. ffTUl) are assumed, together with u; = 10 and t] = T = 10~^. 
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